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Abstract 

We study the collective escape dynamics of a chain of coupled, weakly 
damped nonlinear oscillators from a metastable state over a barrier 
when driven by a thermal heat bath in combination with a weak, glob- 
ally acting periodic perturbation. Optimal parameter choices are iden- 
tified that lead to a drastic enhancement of escape rates as compared 
to a pure noise-assisted situation. We elucidate the speed-up of escape 
in the driven Langevin dynamics by showing that the time-periodic 
external field in combination with the thermal fluctuations triggers an 
instability mechanism of the stationary homogeneous lattice state of 
the system. Perturbations of the latter provided by incoherent ther- 
mal fluctuations grow because of a parametric resonance, leading to 
the formation of spatially localized modes (LMs). Remarkably, the 
LMs persist in spite of continuously impacting thermal noise. The 
average escape time assumes a distinct minimum by either tuning the 
coupling strength and/or the driving frequency. This weak ac-driven 
assisted escape in turn implies a giant speed of the activation rate of 
such thermally driven coupled nonlinear oscillator chains. 

Ever since the seminal work by Kramers (for a comprehensive review see 
Ref. [1]) we witness a continual interest in the dynamics of escape processes of 
single particles, of coupled degrees of freedom or of chains of coupled objects 
out of metastable states. To accomplish the escape the considered objects 
must cross an energetic barrier, separating the local potential minimum from 
a neighboring attracting domain. From the perspective of statistical physics 
mainly the thermally activated escape, based on the permanent interaction of 
the considered system with a heat bath, has been studied p]. The coupling 
to the heat bath causes dissipation and local energy fluctuations and the 
escape process is conditioned on the creation of a rare, optimal fluctuation 
which in turn triggers an escape. To put it differently, an optimal fluctuation 
transfers sufficient energy to the system so that the system is able to statist!- 



cally surmount the energetic bottleneck associated with the transition state. 
Characteristic time-scales of such a process are determined by the inverse of 
corresponding rates of escape out of the domain of attraction. Within this 
topic, numerous extensions of Kramers escape theory and of first passage 
time problems have been widely investigated [H |2] . Early generalizations to 
multi-dimensional systems date back to the late 1960's [3J. This method is 
by now well established and is commonly put to use in biophysical contexts 
and for great many other applications occurring in physics and chemistry 
and related areas [1]-[I1]. 

In order that the system comprised of coupled units may pass through 
a transition state an activation energy Eact has to be concentrated in the 
corresponding critical localized mode (LM). In view of controlling the process 
of barrier crossing we intend to demonstrate that the formation of the critical 
LM can be distinctly accelerated via the application of a weak external ac- 
driving. By use of optimally oscillating barrier configurations it is feasible 
that a far faster escape can be promoted, leading to a drastic enhancement 
of the escape dynamics. Particularly at low temperatures, where the rate of 
thermal barrier crossing is exponentially suppressed, such a scenario can be 
very beneficial. 

Prior studies mainly dealt with the appearance of LMs in damped, driven 
deterministic nonlinear lattice systems [15]- [19]. Furthermore, the sponta- 
neous formation of LMs (breathers) from thermal fluctuations in lattice sys- 
tems, when thermalized with the Nose-method [20] has been demonstrated in 



j . Here we explain LM formation in a stochastic system involving dis- 
sipation in presence of enduring spatio-temporal random forcing. In addition 
a weak external ac-field is applied rendering coherently oscillating barriers. 
The energy is introduced in the lattice coherently in the form of a plane wave 
excitation as the response to the external ac-field and non-coherently through 
thermal fluctuations. We shall demonstrate that the stochastic source and 
the external ac-field conspire to produce such an instability mechanism of 
the stationary fiat state (plane wave) solution yielding a spatially localized 
system state. Most importantly the formed LMs prove to be robust despite 
the continuously impacting thermal forces. 

It should be noted that thermally activated escape of ensembles of non- 
interacting (individual) particles over a metastable potential landscape that 
is additionally subjected to either stochastic or coherent perturbations in the 
form of fluctuations or periodic driving has been studied in the prior litera- 
ture, e.g. see in Ref. [231, [2^ [25] . For a comprehensive overview we refer the 
reader to Ref. [26]. In particular, a resonant activation is observed, i.e., the 
mean escape time (or the rate of escape [27]) attains a minimum (maximum) 
as a function of the correlation time of the fluctuations or the temporal driv- 
ing period of the underlying potential variations. Moreover, the kink drift 
motion induced by oscillating external fields needs to be mentioned in this 
context [28j. Concerning a system of coupled elements the kink-antikink 
nucleation within a (p^ chain model subjected to a deterministic periodic sig- 
nal and uncorrelated noise has been studied in [5]. For optimal noise and 



coupling strength spatiotemporal (array enhanced) stochastic resonance is 
observed in the array of overdamped coupled elements. With the present 
study we focus on the collective nature of the ac-driven escape process of 
interacting weakly damped particles. 

In detail, we study a one-dimensional lattice of damped nonlinear and 
ac-driven coupled oscillators which are subjected additionally to a heat bath 
at temperature T. Throughout the following we shall work with dimension- 
less parameters, as obtained after appropriate scaling of the corresponding 
physical quantities. The coordinate q of each individual nonlinear oscillator 
with a unit mass evolves in a cubic, single well on-site potential of the form 

U{q) = ^q' - Iq'. (1) 

This potential possesses a metastable equilibrium at g^m = 0, corresponding 
to the rest energy Emin = and exhibits a maximum that is located at 
qmax = ^o/a with energy Emax = A_E = a;o/(6a^). Thus, in order for 
particles to escape from the potential well of depth AE over the energy 
barrier and subsequently into the range q > qmax, a sufficient amount of 
energy need to be supplied. The lattice dynamics is governed by the following 
system of coupled Langevin equations 

(in + iqn + ^Iqn - aql + in{t) 

- K [qn+i + qn-i - 2qn] - f sm{ujt + Oq) = . (2) 

The coordinates g„(t) quantify the displacement of the oscillator in the local 
on-site potential U at lattice site n G [1,A^]. The oscillators, referred to as 



"units", are coupled bi-linearly to their neighbors with interaction strength 
n. The friction strength is measured by the parameter 7 and S,n{t) denotes 
a Gaussian distributed thermal, white noise of vanishing mean {C,n{t)) = 0, 
obeying the well-known fluctuation-dissipation relation 

with /cb denoting the Boltzmann constant. A homogeneous external periodic 
modulation field of amplitude /, frequency u and phase 6*0 globally acts upon 
the system. In this work we use periodic boundary conditions according 
to qN+i = Qi and fix the parameters of the potential as follows: uJq = 2 
and a = 1, yielding AE = 4/3. A deterministic escape scenario in the 
conservative, undriven limit of system ([2]) has been explored by us in [2^150] . 
To analyze the nonlinear character of the solutions of Eq. ([2]) we first 
discard the noise (^„ = 0) and derive a nonlinear damped and driven dis- 
crete Schrodinger equation for the slowly varying envelope solution, Un{t), 
following the reasoning in [31], i.e.. 



2iuJo iin + i^UJoUn + K. [Un+1 + M„_i - 2m„] 

+ a\un\'^Un + -f exp[-iAujt + 6o] = 0, (4) 

with the nonlinearity parameter reading a = 10a^/3u;g and Au = u — uq. 
For the amplitude uq of a spatially homogeneous solution of Eq. (jlj of the 
form 

Unit) = uq exp[— 'j(Ati;t + 9o)] + c.c. (5) 
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Figure 1: (color online) Response of the amplitude Uq with respect to the 
frequency u of the driving force for two fixed values of the driving amplitude 
as indicated in the plot. The damping constant is 7 = 0.1. The inset shows 
the response of the amplitude uq with respect to the driving amplitude / for 
a fixed frequency u = 1.295. The remaining parameter values are Uq = 2 
and a = 1. 



one obtains 



ul 



;f- 



(6) 



{2ujqAuj + aul) + j'^ujq 

In Fig. [T] we depict the amplitude Uq of the response versus the driving 
frequency curve for two different values of the driving amplitude /. At a 
bifurcation point a "jump" resonance related with a saddle-node bifurcation 
occurs and in certain range of the driving frequency multistability exists. In 
comparison for the larger driving amplitude, / = 0.2, the bifurcation point 
for the "jump" resonance occurs at a lower frequency value than for the 
driving with / = 0.15. Moreover, in the former case the system responds 
overall with higher amplitudes Uq than in the latter. Notice that the system 



responds with large amplitude only within a frequency window and large 
amplitudes are obtained for a driving frequency lying below the band of 
linear frequencies, viz. to < ujq = ^J2 = 1.414.... Similarly, for the response 
of the amplitude with regard to the driving strength / multistable solutions 
are possible as depicted with the inset in Fig. [H In order to investigate the 
stability of the homogeneous solution of Eq. ([2]) we use 



qn{t;k = 0) = x{t) = Uoe-'^^'+''>^ 

ul + c.c. . (7) 



a 



^0 L 



9 _ l^-2i(a;t+eo) 

3^^ 



and write with respect to the spatial perturbations A„: qn{t) = x{t) + A„(t). 
Since we impose periodic boundary conditions the Fourier-series expansion 
^n(t) = J2k exp(zfcra)sfc(t) can be used to yield an equation for the mode 
amplitudes Sk, i.e., 

Sk + ISk + iolsk - 4aMo cos(u;t + 6o)sk = , (8) 

where we discarded a higher harmonics and introduced ul = Uq+Ak sin^ {k/2)- 
8{auQ/uJo) . Setting r = ujt/2 and Sjt(t) = ffc(t) exp(— 7t/2) one derives a 
Mathieu equation 

Vk + [A- 2Q cos(2t + 2^o)]f fc = , (9) 

with the parameters A = {2uJk/uj)'^ — {'y/uj)'^ and Q = Sauo/u"^. If it holds that 
y/A ~ /, with / denoting a positive integer number, the Mathieu equation 
allows for parametric resonance [32], [33]. The extension of the resonance 
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Figure 2: (color online) Instability bands for different coupling strengths 
K = 2, 1, 0.5, 0.3, 0.2 (decreasing from left to right) are displayed. The relation 
between the amplitude of the external ac-field / and the amplitude of the 
homogeneous solution uq is the one displayed with the inset in Fig. [TJ The 
horizontal dashed line at / = 0.15 intersects each instability band very close 
to its bottom, the position of which determines the respective critical wave 
number kc (see also text). 



regions is determined by the ratio Q/A; for the primary resonance, A ~ 1, it 
is given by 

iA-lf<Q\ (10) 

For the parameter set corresponding to the line shown in the inset in Fig. [1] 
(determining the relationship between the amplitude / of the external ac- 
field and the amplitude of the homogeneous solution uq) the instability bands 
for different values of the coupling strength k are depicted in Fig. [2l For the 
onset of parametric resonance the driving amplitude / has to exceed the value 
of the bifurcation point, i.e. fc ^ 0.1408 related with the "jump" resonance. 



regardless of the value of k. The position of the bottom of the instability 
band, determining the critical unstable wave number kc, shifts towards lower 
fc— values with increasing coupling strength k. For a chosen field strength 
/ = 0.15, that lies just above /c, one expects that the LMs of distinct wave 
length, determined by Ac = 271 /k^ become excited (cf. Fig. [2]). We infer 
from Fig. [2] that the wave length of a LM increases with increasing coupling 
K. This is verified in Fig. [3] showing the spatio-temporal evolution of the 
amplitudes qnif) for couplings k = 0.5 and k = 2. The Langevin equations 
were numerically integrated using a two-step Heun stochastic solver. In all 
our simulations the initial chain configuration is represented by gn(0) = a;(0) 
and Pn(0) = with the homogeneous solution (plane wave) x(0) given in 
(ITj). We note the formation of a LM of certain wave length arising from the 
homogeneous state after a short time span (after t ~ 60) and we note that 
the period duration for oscillations near the bottom of the potential is around 
2-n/ujQ ~ 4.4. 

Regarding the energy relation within the stationary fiat state (where each 
unit contains the same amount of initial energy) we monitored the temporal 
evolution of the energy of one unit 

En = \vl + f/(gn) , (11) 

and the corresponding field energy 

E field = -f sm{ujt + 9o)qn , (12) 

without coupling the chain to the heat bath for a force amplitude / = 0.15. 
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Figure 3: Spatio-temporal pattern of the solutions qn{t) for a lattice consist- 
ing of A^ = 100 sites and chosen coupling strength at k = 0.5 (top panel), 
K, = 2 (central panel) for the same realization of Gaussian white noise with 
thermal energy ksT = 0.001 x AE. Bottom panel: Same as in the central 
panel but now for a 50-times larger thermal energy: ksT = 0.05 x AE. The 
remaining parameter values are set at / = 0.15, u = 1.295, ^o = and 
friction 7 = 0.1. 



In this stationary case the field energy performs small-amphtude oscillations 
around a mean value of Efieu = 0.04 = 0.03 x /S.E = 0.03 x 4/3, while the 
mean of the energy of one unit is En = 0.78 = 0.585 x AE (not shown). Thus 
the gain of energy, determined by the ratio En/Efi^iii, amounts to a remark- 
able high value of 19.5. To retain this relation upon lowering (increasing) the 
damping 7 a lower (higher) driving strength / is necessary while the "jump" 
resonance frequency attains a lower (higher) value according to Eq. ([6]). 

The stochastic term provides perturbations of all wave numbers and a 
pattern emerges from the homogeneous fiat state. That is, due to the effect of 
parametric resonance perturbations provided by the thermal noise grow and 
induce a LM consisting of several humps. The fastest growing perturbations 
are those associated with the critical wave number kc (see also [34J). Each 
of these humps resembles the hairpin shape of the transition state as the 
critical escape configuration possessing an energy Eact through which the 
coupled units have to pass in order to cross the barrier [30]. The robustness 
of the LMs is remarkable: a LM is sustained despite continuously impacting 
thermal noise of strengths up to values ksT < 0.2 x AE. Moreover, the 
formed pattern maintain their distinct wave length Ac = 211 /kc (see Fig. |3]). 

We note that upon increasing the noise strength the growth rate of the 
humps becomes enhanced, being refiected in the statistics of the barrier cross- 
ing of the chain in the presence of weak ac-driving. The amplitude and fre- 
quency of the latter are chosen such that the dynamics exhibits parametric 
resonance. The dependence of the mean escape time of the chain on the in- 
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Figure 4: The mean escape time of the chain versus the mean injected energy 
E = E field + ksT measured in units of AE with fixed field energy Eji^id = 
0.03 X AE provided by an external modulation field with u = 1.295, ^o = 
and / = 0.15. Here we vary the thermal energy E thermal = ksT. The inset 
depicts the unforced case with / = 0. The remaining parameter values are 
N = 100, K = 0.28 and 7 = 0.1. 



jected average energy E = Efi^d + Ethermau with Ethermai = ^bT (measured 
in units of the barrier energy AE) is displayed in Fig. |H The thermal energy 
Ethermai 1 Supplied nou-cohereutly by the heat bath, is varied Wiihm the range 
[(10-^-0.11) X AE]. 

The average of the escape times was performed over 500 realizations of the 
thermal noise. In this context the random escape time of a unit is defined as 
the time instant when the unit passes through the value g = 20 far beyond the 
potential barrier. Thus, no likely recrossing back into the potential valley can 
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occur [221 EOj- The escape time of the chain is then determined by the average 
of the escape times of its units. We notice that the underlying irregular 
dynamics serves for self- averaging and thus the choice of the phase of the 
coherent, external forcing, 6q, does not affect the mean escape time. In the 
forced as well as unforced case there occurs a rather rapid decay of Tesc with 
growing Ethermai = ^bT at low temperatures. This effect weakens gradually 
upon further increasing ksT. Most strikingly, for the forced system the 
escape times become drastically shortened in comparison with the unforced 
case with / = 0. Moreover, for the forced system escape takes place also at 
very low temperatures for which in the undriven case not even the escape 
of a single unit has been observed during the simulation time (taken here as 
t = 10^) implying a giant enhancement of the rate of escape as compared to 
the purely thermal noise driven rate. 

Upon exploring the optimal escape route we investigated the influence 
of the coupling strength k on the average escape time. Our numerical find- 
ings are summarized in Fig. [51 The mean escape time exhibits a resonance 
structure, viz. there exists an optimal coupling strength [Kres — 0.28) for 
which the escape assumes a minimum. Upon lowering k, < n^es we notice 
a drastic rise of the escape time while for k, > Hres the graph exhibits only 
a moderately growing slope with growing coupling strength k. We empha- 
size the collective nature of this resonance effect which here occurs for finite 
interaction strength k 7^ 0. In the limit k, —>■ the mean escape time of non- 
interacting, individual particles assumes for this parameter set an extreme 
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Figure 5: The mean escape time versus the couphng strength k exhibits a 
resonance structure. The inset displays the ratio R, defined in Eq. flT3|) . as 
a function of n. The remaining parameter values are given by A^ = 100, 
ksT = 0.05 X AE, f = 0.15, u = 1.295, ^q = and 7 = 0.1. 



large value, implying a vanishingly small escape rate. 

To explain the occurrence of the resonance structure in Fig. [5] we recall 
that the wave length, Ac = 2iT/kc, of the arising LMs on the lattice is deter- 
mined by the critical wave number kc = kc{K,, f) (cf. Fig. [2]). The number of 
humps contained in a LM, Nh-, can be attributed to kc as: XcNh = 2'K/kcNh = 
N. The number of humps (besides their height and width) regulates how the 
mean energy injected via the coherent external field and the incoherent ther- 
mal noise is shared among them. Supposing that the whole lattice can be 
divided into an array of segments, where each of them supports a single local- 
ized hump, the energy of one segment is given by Eg = E/Nh = 27TE/{kcN). 
Appropriate conditions for successful escape are provided when the energy 
contained in each segment, E^, is close to the activation energy, Eact of the 
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Figure 6: The mean escape time as a function of the driving frequency. The 
thermal energy is ksT = 0.05 x AE and the driving amphtude is / = 0.15. 
The remaining parameter values are A^ = 100, uj = 1.295, k = 0.3, / = 0.15 
and 7 = 0.1. For comparison we note that Langer's theory [3] yields for the 
mean escape time in the unforced case, / = 0, the value T^sc — 553 x 10^ 
underpinning the drastic speed-up effect in our forced case. 



critical escape configuration [30] . The efficiency of energy localization is then 
determined by the ratio 

R = Es/Eact- (13) 

The activation energy as a function of the coupling strength satisfies (we 
recall that we use a dimensionless formulation) the relation Eact = (1 + 3.54 x 
K,)AE [5U]. Keeping the injected energy E fixed and given value of k^ we 
obtain R. In the inset of Fig. [5] the ratio R is plotted as a function of the 
coupling strength k. The plot indeed exhibits a maximum at k = 0.28, which 
conffims the finding of the resonance found for the mean escape time versus 
coupling strength as depicted in Fig. [51 Concerning the critical localized 
mode through which a lattice state has to pass through in order to escape 
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over the potential barrier we remark that for comparatively low coupling 
strengths (k < 0.6) the effect of discreteness in the lattice system is still such 
pronounced that this critical localized mode is indeed represented by a thin 
hairpin-shaped configuration involving effectively one lattice unit of large 
amplitude to either side of which the amplitude pattern decays extremely 
rapidly (for more details see [3D])- 

Next we study the role of the angular driving frequency of the external 
modulation field, see in Fig. [61 The escape time as a function of the angular 
frequency likewise exhibits a resonance structure and there exists an optimal 
frequency for which the average escape time assumes a minimum. This is 
reminiscent of the phenomenon of resonant activation found for the thermally 
activated escape of noninteracting particles surmounting oscillating barriers 
[23l [211 [251 ISSl [27]. In our case of a nonlinear chain composed of coupled 
units, however, this 'resonant activation' within a frequency window nicely 
correlates with the systems' gain of energy that is coherently supplied by the 
applied ac-field in this very same frequency interval (note the corresponding 
frequency window in the nonlinear frequency response graph associated with 
large amplitudes in Fig. [T|). Therefore, tuning the frequency u; at a fixed 
interaction strength k allows to optimize further the mean escape time. The 
minimal escape scenario thus requires an optimal tuning both in coupling 
strength and ac-driving frequency uj. 

In summary, we have presented a drastic speed-up mechanism of the 
thermal noise driven barrier crossing of a coupled damped nonlinear oscil- 
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lator chain under the impact of a weak external ac-field. With appropriate 
parameter values of the latter and in the presence of thermal noise an in- 
stability mechanism is initiated due to which LMs arise from stationary flat 
state solutions in the lattice dynamics. Humps of the LMs are rapidly driven 
through the transition state thus accelerating the escape over the situation 
with purely thermally assisted escape. Interestingly, the LMs are sustained 
up to fairly high noise levels corresponding to ksT ~ 0.2 x AE. 

The findings of our study can be applied for the control of the rate of bar- 
rier crossing of oscillator chains. With such chains providing the archetype 
model for nonlinear collective transport of matter, charge and energy in abun- 
dant low-dimensional systems in physics, biology and chemistry this speed-up 
scenario of thermally driven collective escape over potential barriers might 
well be put to constructive use in a variety of potential applications. 
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